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1 Problem 

Deduce the velocity distribution of steady flow of an incompressible fluid of density p and 
viscosity rj between two parallel, coaxial annular plates of inner radii 7*1, outer radii r2 and 
separation h when pressure difference AP is applied between the inner and outer radii. 

As an exact solution of the Navier-Stokes equation appears to be difficult, it suffices to 
give an approximate solution assuming that the velocity is purely radial, v = v(r, z)r, in 
a cylindrical coordinate system (r, <p, z) whose z axis coincides with that of the two annuli. 
Deduce a condition for validity of the approximation. 

This problem arises, for example, in considerations of a rotary joint between two sections 
of a pipe. Here, we ignore the extra complication of the effect of the rotation of one of the 
annuli on the fluid flow. 



2 Solution 

For an incompressible fluid, the velocity distribution obeys the continuity equation 

V-v = 0, (f) 

in which case the Navier-Stokes equation for steady, viscous flow is 

p(v • V)v = -VP + r/V 2 v. (2) 

There are only three examples in which analytic solutions to this equation have been obtained 
when the nonlinear term (v ■ V)v is nonvanishing [T[|. 

We first review the simpler case of two-dimensional flow between parallel plates in sec. 2.f , 
and then take up the case of radial flow in sec. 2.2. We will find an analytic solution to the 
nonlinear Navier-Stokes equation (0) for radial flow, but this solution cannot satisfy the the 
boundary conditions, 

v(z = 0) = = v(z = h), (3) 

that the flow velocity vanish next to the plates. However, in the linear approximation to 
eq. (HD we obtain an analytic form for the radial flow between two annular plates. 



2.1 Two-Dimensional Flow between Parallel Plates 

For guidance, we recall that an analytic solution is readily obtained for the related problem 
of two-dimensional viscous flow between two parallel plates. For example, suppose that the 
plates are at the planes z = and z = h, and that the flow is in the x direction, i.e., 
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v = v(x,z)x.. The equation of continuity ([[]) then tells us that dv/dx = 0, so that the 
velocity is a function of z only, 

v = v(z)'k. (4) 

The z component of the Navier-Stokes equation (g) reduces to dP/dz = 0, so that the 
pressure is a function of x only. The x component of eq. (0) is 

dP{x) _ d 2 v(z) 
dx ~ V dz 2 ■ (5) 

Since the lefthand side is a function of x, and the righthand side is a function of z, equation 
(H) can be satisfied only if both sides are constant. Supposing that the pressure decreases 
with increasing x, we write 

dP AP 
— — = constant = — — > 0. (6) 
ox Ax 

Using the boundary conditions (|3|), we quickly find that 

. . APzih-z) z ( z\ 

where the average velocity v is given by 

1 r h AP h 2 

v(z)dz 



h Jo Ax 12r) 

2.2 Radial Flow between Parallel Annular Plates 

Returning to the problem of radial flow between two annular plates, we seek a solution 
in which the velocity is purely radial, v = v(r,z)r. The continuity equation (|l|) for this 
hypothesis tells us that 

1 d(rv) , , 

~^ = 0, 9 
r or 

so that 

v = M*. (10) 
r 

Following the example of two-dimensional flow between parallel plates, we expect a parabolic 
profile in z as in eq. (0), 

f(z)<xz(h-z), (11) 

which satisfies the boundary conditions (|3|). 

Using the trial solution (|T0|) , the z component of the Navier-Stokes equation (0) again 
tells us that the pressure must be independent of z: P = P(r). The radial component of 
eq. (H) yields the nonlinear form 

W 2( % + P f = r id -f. (12) 
az z ar 
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The hoped-for separation of this equation can only be achieved if f(z) = F is constant, 
which requires the pressure profile to be P(r) = A — pF 2 /2r 2 . The boundary conditions (||) 
cannot be satisfied by this solution. Further, this solution exists only for the case that the 
pressure is increasing with increasing radius. The fluid flow must be then be inward, so the 
constant F must be negative. The Navier-Stokes equation is not time-reversal invariant due 
to the dissipation of energy associated with the viscosity, and so reversing the velocity of a 
solution does not, in general, lead to another solution. 

While we have obtained an analytic solution to the nonlinear Navier-Stokes equation (fj), 
it is not a solution to the problem of radial flow between two annuli. It is hard to imagine a 
physical problem involving steady, radially inward flow of a long tube of fluid, to which the 
solution could apply. 

Instead of an exact solution, we are led to seek an approximate solution in which the 
nonlinear term f 2 of eq. fll~2| ) can be ignored. In this case, the differential equation takes the 
separable form 

d 2 f rdP 

-— = — — = constant. (13) 
az A T) dr 

Following eq. (|7|) we write the solution for / that satisfies the boundary conditions @ as 

^^K 1 -!)- (14 > 

where / is the average of f(z) over the interval < z < h. The part of eq. ( |T3"D that describes 
the pressure leads to the solution 

Pi In r 2 /r + P 2 In r/n 

P(r) = hT^ ' (15) 



where Pi = Piji). Plugging the solutions (|i~4l) and ( |l5j) back into eq. fll3l) , we find that 



j_ h 2 AP 
12r) In r 2 /r\ ' 

where AP = P\ — P 2 . Hence, the flow velocity is 

. . z{h-z)AP„ 
v(r,z)= y ' r, 17 
lr\r In r 2 /ri 

whose average with respect to z is v(r) = f/r. As with all solutions to the linearized 
Navier-Stokes equation, the velocity is independent of the density. 

For the approximate solution ( |T7D to be valid, the term f 2 w f 2 must be small in eq. (|T2|), 
which requires 

ph 4 AP . . 

144r] 2 r 2 lnr2/ri 

When this condition is not satisfied, the solution must include velocity components in the 
z direction that are significant near the inner and outer radii, while the flow pattern at 
intermediate radii could be reasonably well described by eq. fllTD . 

If one of the annuli is rotating at angular velocity u, the radial flow velocity should still 
be given approximately by eq. (0) so long as ujr 2 ^ v{r 2 ) = f /r 2 . 
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